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Abstract 

We construct the effective chiral Lagrangian involving the D-mesons and Goldstone bosons 
at next-to-leading order taking into account strong as well as electromagnetic interactions. 
This allows us to disentangle — to leading order in isospin violation — the electromagnetic 
and the strong contribution to the Z?-meson mass differences. In addition, we also apply 
the interaction to the decay £>* (2317) — > D s ir° under the assumption that the Z3* (2317) 
is a hadronic molecule. We find (180 ± 110) keV for the decay width r(D* (2317) -> D s n°) 
- consistent with currently existing experimental constraints as well as previous theoretical 
investigations. The result provides further evidence that this decay width can serve as a criterion 
for testing the nature of the D* (2317). 

PACS: 12.39.Fe; 12.39.Mk; 13.25.Ft; 14.40.Lb 
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1 Introduction 

In recent years, many heavy mesons with open or hidden charm were discovered, which contribute 
to the revival of hadron spectroscopy (for recent reviews, see [1]). One outstanding example among 
them is the D* (2317) discovered by the BABAR Collaboration in the D s ir final state [2]. The 
measured mass of the Z?* (2317) which is 2317.8 ± 0.6 MeV [3] is much lower than that predicted 
in many quark models. An appealing alternative is that the D* (2317) is a hadronic molecule, 
which means that it owes its existence to meson-meson dynamics [4, 5]. In this work we exploit 
further this idea. Note, in Refs. [6] it was argued that a molecular interpretation of the D* (2317) 
(and its vector counter part) is at variance with heavy quark effective field theory. However, this 
conclusion is based on the assumption that the decay of a hadronic molecule is proportional to the 
molecular wave function at the origin — in Ref. [7] it is shown that this assumption is not justified 
for the decay of hadronic molecules. 

The mass alone is not a signal for a molecule, as stressed, e.g., in Ref. [8]. A consistent treatment 
of the mass and various decays is required — see also Ref. [9]. The difference from the molecular 
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state will presumably be revealed in the decay pattern into various channels. However, no branching 
ratio of the D* (2317) has been reported accurately. The only experimental constraints are upper 
limits for the ratios of some other decay channels to the Z)* (2317) + — > Dfir . For instance [3], 

r(Z? s * (23i7)+ ^(2112)+ 7 ) 

r( J D s * (2317)+ - D+ttO) 1 j 

The D* (2317) can be dynamically generated from Goldstone boson-D-meson scattering as a 
hadronic molecule using unitarized amplitudes from chiral pertrubation theory [10, 11, 12, 14]. 
The width of the isospin violating decay D* (2317) + — ► D^tt was estimated to be about 8.7 keV 
in Ref. [12] by considering only the Tr°-r] mixing in the final state. However, Faessler et al. pointed 
out that the mass differences between neutral and charged kaons and L>-mesons give an important 
contribution [13], which was confirmed later in Ref. [14]. In Refs. [11, 14] also subleading operators 
were studied. Note that isospin symmetry violation in hadronic physics has two different sources: 
one originates from the mass difference of the light u and d quarks, and the other one stems 
from the electromagnetic (e.m.) interaction. While the TT°-r] mixing and part of the meson mass 
differences account for the former one, the effect of the latter on the decay -D* (2317) + — > -D+7T 
will be investigated here for the first time. 

The first step is to construct the interaction Lagrangian to next-to-leading order in the chiral 
expansion. Based on this we can disentangle the e.m. and the strong contribution to the D-meson 
mass difference. As will be demonstrated, the Lagrangian also links these mass differences directly 
to the isospin-violating Goldstone boson-D-meson scattering amplitudes, in full analogy to the 
case of ttN scattering and the proton-neutron mass difference [15, 16]. We can therefore calculate 
the decay width of the D* (2317) + — > D+ir° — within the molecular picture utilizing the chiral 
Lagrangian up to the next-to-leading order (NLO) 0(p 2 ), where p denotes a small parameter - 
with a small number of free parameters. Note, in this paper for the first time both the strong and 
the e.m. contributions to the decay are incorporated systematically. 



2 Lagrangians at next-to-leading order 

The scattering between the Goldstone bosons and .D-mesons is similar to the case for pion-nucleon 
scattering (for reviews, see Refs. [17, 18]) because the D-mesons have heavy masses which do not 
vanish in the chiral limit* 1 . We count the D-meson masses (~ 1.9 GeV) as order 0(A X ) ~ O(p ) 
where A x ~ 1 GeV. Hence the leading order terms in the chiral Lagrangian are of O(p), and the 
NLO terms are of 0(p 2 ). 

The leading order Lagrangian is just the kinetic energy term of the heavy mesons [19] 

= V^DV^rf - m 2 D DD^ (2) 

with D = (D°,D + ,Df) denoting the D-mesons, and the covariant derivative being 

v„ = ^ + r M , 

r ^ = o ( ujd v u + ud^) ' ( 3 ) 

where 




U = exp -V - , u z = U. (4) 



#1 In this work we consider the SU(3) chiral limit, m u , md, m s — > 0. 
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The Goldstone boson fields are collected in the matrix 



V2" 1 V6 

K- K° 



(5) 



We now consider the NLO chiral Lagrangian describing the interactions of the pseudoscalar 
charm mesons with the Goldstone bosons. Considering the heavy mesons as matter fields, similarly 
to the pion-nucleon sector [16], the strong part is 

££l = D{-h (x + )-h 1 x + + h 2 (u II u»)-h 3 u fl u»)D 

+V ll D{h 4 {u»u v ) - h 5 {u»,u u } - h & [u»,u v ])V v D. (6) 

The chiral symmetry breaking terms, i.e. ho and h± terms, have been introduced before [20, 11]. 
The /12 and /13 terms were introduced in Ref. [11]. We stress that the contributions of the /15 
and /13 terms to s-wave amplitudes differ only to order 0{p/mri). However, we still keep them 
in our covariant formalism for in this way we have an additional tool to estimate the theoretical 
uncertainty — see section 4.2. The electromagnetic part is 

41 = FlD [<7 (Ql - Q 2 _) + 91 (Ql -Ql)+ g 2 Q+ (Q+) + 33 (Q+) 2 } D , (7) 

where 

X+ = u'xu' + uxu, 
X+ = X+ ~ 3 <X+) 1 

Q± = \ (u*Qu ± uQu^ . (8) 

The quark mass matrix and the D-meson charge matrix are diagonal 

X = IB ■ di&g {m u ,m d ,m s } , Q = e ■ diag {0, 1, 1} , (9) 

in terms of B = \ (Q\qq\Q)\/ F% and the elementary charge e. Further, F w is the pion decay constant. 
The unknown coefficients hi (i = 0, 1, 6) and Qi (i = 0, 1, 2, 3) in Eqs. (6,7) are the so-called low 
energy constants (LECs). As we will show in the next section, h\ and a linear combination of g$ 
and g 2 , namely go + 2g 2 , can be determined from the mass differences among the .D-mesons. Since 
(Q+) = 2e, the 53 term in C e .-L. only gives an overall e.m. mass shift of the D-mesons, and hence 
can be absorbed in the bare masses. The g\ term contains only isospin symmetric e.m. interaction 
and is irrelevant here. Terms with one more flavor trace in the strong interaction Lagrangian are 
suppressed in the large Nq limit of QCD [21]. We therefore follow Ref. [14] and drop the ho, h 2 and 
h^. Formally the h& term is of 0(p 2 ). However, due to the commutator structure, it is suppressed 
by one order, see Appendix A. We are therefore left with only two free, active parameters, both 
isospin conserving, namely /13 and h§ . We will investigate their effect on the isospin violating decay 
of the D* below. 

One should note that the 7r° and 77 in Eq. (5) are not mass eigenstates because of the 7r°-r/ 
mixing. The mass eigenstates are defined as 

7T° = 7r°cose 7r o J) + ?7sine 7r o^, 
fj = — 7T° sin e n o v + r/ cos e 7T o v , (10) 
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where e n o v is the well-known 7r°-r] mixing angle, which reads to leading order 

V3m d -m u 

Mj? = — ; — I 11 ) 

4 m s — m 

with m = (m u + md)/2 the average mass of the u and d quarks. 



3 D— meson mass differences 

The terms which can contribute to the mass differences among the D + , D° and Df mesons come 
from the Lagrangian of order 0(p 2 ), Eqs. (6,7). Only three among all the terms contribute. We 
find 



m D0 - m D+ = hX + g 

2 



m 2 D+ -m 2 D+ = /»(!-£), (12) 



where we define 



A = — ~ = — e *°v ( 14 ) 



h = ABh^m.-m), g = F 2 e 2 (g + 2g 2 ). (13) 

The strength of isospin violation due to quark mass effects is encoded in the parameter A that is 
connected to e^o^ through 

m d - m u 
m s -rh y/jf 

A recent analysis of p-uj mixing in chiral perturbation theory gives 1/A = 42 ± 4 [22], correspond- 
ingly we have A = 0.024 ± 0.002. Using the masses for the D-mesons [3], m D o = 1864.84 ± 0.17 
MeV, m D + = 1869.62 ± 0.20 MeV, and m n+ = 1968.49 ± 0.34 MeV, we find 

h = (384.1 ± 2.5) x 10 3 MeV 2 , g = (4 ± 1) x 10 3 MeV 2 , (15) 

where the largest uncertainty comes from the masses of the D-mesons. Note that here only the 
uncertainties from the experimental inputs are considered. For a discussion of the theoretical 
uncertainty, see Section 4.2. Then the dimensionless LECs h\ and go + 2g 2 can be determined as 

hi = 0.42 ± 0.00, g + 2g 2 = 11 ±3, (16) 

where we use B(m s — rh) = (M^ + M|. + ) /2 — M 2 . The basic assumption in setting up an 
effective field theory is the naturalness of the low energy constants — especially the dimensionless 
coefficient h\ should be of order one. This is indeed the case for we find h\ = 0.42. A naturalness 
estimate for go and g 2 comes from requiring that the contribution of the corresponding operators 
to the .D-meson mass shift should be of the order of a typical virtual photon loop [23], thus 

e 2 4~(^) 2 < (17) 

This leads to gi ~ 4 as a natural estimate of the order of magnitude, compatible with the value 
determined for go + 2g 2 . 

The parameter h, fixed from the amount of 51/(3) violation encoded in the mass difference 
between the Df and the D + (see Eq. (12)), controlls also the strong part of the D° and D + mass 
difference. Therefore, the electromagnetic contribution to this mass difference, which is given by 
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the g term, can be extracted from data. This is different to the case of, e.g., nucleons, where the 
operator structure is more complicated. We therefore get 

(m D +-m D o) stl = (2.5 ±0.2) MeV, (18) 
(m D+ - m D0 ) e m = (2.3 ± 0.6) MeV, (19) 

where the first equation refers to the strong contribution to the mass difference and the second to 
its e.m. counterpart. Note, contrary to what is common for nucleons as well as kaons, here the 
electromagnetic and the strong effects enter with the same sign. This is a direct consequence of 
the different quark content of the states. 

The strong and e.m. mass differences of rri£)+ — m^o are consistent with those determined long 
time ago by Gasser and Leutwyler using a simple quark model ansatz [24], which are 3.3 ±0.9 MeV 
and 1.7 ± 0.5 MeV, respectively. 



4 Width of the £^(2317) 

For studying isospin violating decays, it is better to work in the particle basis. For the scalar 
charm-strange sector, there are four channels involving a D-meson and a Goldstone boson: D°K + , 
D + K°, D+r) and D+tt . One can expect that the isospin violating contributions to the mass of the 
D*q(2317) is negligible, therefore we only consider them for calculating the isospin violating decay 
width. The strong contribution to this decay is given in terms of the known 7r°~T] mixing angle 
and the meson mass differences. The only non-vanishing e.m. contribution is from the transition 
D°K+ D+ir° 

V d™k + ^dU° = -^(90 + 2g 2 )e 2 . (20) 

Especially, the amplitude for Dfrj —> D^n vanishes. The linear combination of LECs (go ± 2g 2 ) 
has been determined from the D-meson mass differences in the previous section. Thus, all relevant 
isospin violating interactions are fixed from data. 

4.1 Unitarization of the scattering amplitudes at next-to-leading order 

A unitarization procedure was proposed in Ref. [25] which can be used for any finite order in the 
chiral expansion. Similar to the case for pion-nucleon scattering, in our case, up to NLO there 
is no loop contribution. Hence we obtain the following T-matrix equation after matching to the 
chiral expansion at NLO 

T( S ) = V( S )[l-G(s).V(s)}- 1 , (21) 

with V(s) = Vlo(s) + Vnlo(s) the sum of the S-wave scattering amplitudes of the LO and NLO 
orders [25]. G(s) is a diagonal matrix with the diagonal element given by the two-meson loop 
integral [25, 26] 

G(s) u = 



d 4 q 



(2vr) 4 ( g 2 - m \ + ie) [(P - qf -m 2 2 + ie] 

,2 ,^2 ,^2 , - ™2 



1 f , mk m? — mi + s , mi a r . , 9 9 , 
= 1 ^{^)±ln^± ^ ln-| + -[ln (g -^ + mI + .) 

— ln(— s + m\ — m 2 .+ a) + ln(s + m\ — m\ + a) — ln(— s — m\ + m\ + a)] } , (22) 

where a(fi) is the subtraction constant, [i denotes the scale of the dimensional regularization, 

1 /2 

and o = {[s — (mi + rri2) 2 ][s — (m\ — m2) 2 ]} . In our analysis we use the somewhat lengthy 
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LO 7r — 77 mixing Mass differences 
149.4 DTO 69T7 



Table 1: Decay widths of the D* (2317) — > -D s 7r° with LO amplitudes. All units of the decay 
widths are in keV. 

expression given above, for it allows for a straightforward analytic continuation into the complex 
plain, contrary to more compact representations that are applicable in a particular parameter 
space only. 

The subtraction constant a(n) is determined by fitting the mass of the D* (2317) using the LO 
Lagrangian. It turns out to be a(fi = 1 GeV) = —1.846 for reproducing mp* = 2317.8 MeV [3]. 
Probably by accident, it coincides exactly with that obtained from matching the value of the loop 
function at the threshold of the D and K calculated by Eq. (22) with that calculated by using a 
3-momentum cut-off g max = rn p [12]. 

4.2 Results 

In calculations, we take physical values of all the meson masses and the pion decay constant, as 
listed in the following [3]: 

F n = 92.42 MeV, M^o = 134.98 MeV, M 7r + = 139.57 MeV, 

M K o = 497.65 ± 0.02 MeV, M K+ = 493.68 ± 0.02 MeV, 

m D o = 1864.84 ± 0.17 MeV, m D + = 1869.62 ± 0.20 MeV, 

m D + = 1968.49 ± 0.34 MeV, M r) = 547.51 ± 0.18 MeV. (23) 

After unitarization, the hadronic molecule D* (2317) appears as a pole in the second Riemann 
sheet at {tud* — ^T(Z)* — > D s 7r°)/2) . Denoting the three-momentum of one particle in the center- 
of-mass frame of channel i by fcj, the second Riemann sheet is specified by Im k D +^ < and Im fcj > 
(i = D°K + , D + K°, Dfrj). First, let us focus on the results considering the LO amplitudes only. 
There are two different kinds of contributions to the isospin violating decay width at leading order. 
One is from the 7r°-ry mixing and the other one is from the mass differences between charged and 
neutral kaons and D-mesons, which predominantly enters through an isospin violating contribution 
to the loop function of Eq. (22). However, e n o v = l/y/3B(m,d — m u )/ (M^ — M^ ) is suppressed 
by M^/M 2 — in SU(2) chiral perturbation theory, where the strange quark is also viewed as 
heavy, this operator appears only at NNLO, even an order below those given in Eqs. (6) and (7). 
Therefore, one can expect that the mass differences give a larger contribution. The results confirm 
this expectation as shown in the second and third column in Table 1, corresponding to the widths 
considering only the 7T°— 77 mixing and meson mass differences, respectively. Furthermore, similar 
to Refs. [13, 14], in our calculation the interference between these two kinds of contributions are 
constructive, giving rise to a width of about 150 keV — see the first column of Table 1. However, 
when it comes to a quantitative comparison, the result for the width to leading order of Ref. [14] is 
smaller by a factor of two — a direct comparison with the more phenomenological work of Ref. [13] 
is not possible. The difference can be traced to differences in the input parameters and a different 
method to fix the subtraction constant a(fi) of Eq. (22). Those differences should be of higher 
order. Thus the spread in the reported results calls for a calculation to next-to-leading order in 
the chiral expansion, c.f. Ref. [14], together with an analysis of the uncertainties. 

We take [—1, 1] as a natural range for the dimensionless parameter h' 5 = h^/m 2 D0 — note that 
for h' 5 = ±1 the contribution of the /15 term to the D°K + — > D°K + scattering amplitude is of the 
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same order as the leading one. The subtraction constant is kept fixed to a(l GeV) = —1.846. The 
value of ti3 is then determined from fitting the pole position in the second Riemann sheet to the 
mass of the scalar charm meson mD* s0 = 2317.8 ±0.6 MeV in the isospin symmetric case. For each 
value of h' 5 , there is a corresponding /13, as shown in Fig. 1. 

For each value of h' 5 the resulting widths are plotted in Fig. 2. In the figure, the solid line 
is the result for producing the mass of the -D* (2317) at 2317.8 MeV with both the strong and 
e.m. isospin violating contributions and using central values for all the parameters. The dashed 
curve represents the result without e.m. contributions. The uncertainties from the experimental 
inputs are reflected in the shaded band. To estimate the theoretical uncertainty observe that, 
as stressed above, the effect of the /13 and the /15 term differ at O{p/mo). Therefore, since one 
of the two is fixed already from the mass of the D s q, the dependence of the width on /15 is a 
measure of (some) NNLO effects. We may read off the figure directly a spread of about 50 keV 
around the central value of 180 keV induced by the variation of /15. To be on the safe side we 
take as the theoretical uncertainty of our calculation twice this spread. Another way to estimate 
the theoretical uncertainty is to use 2(Mx/A x ) 2 ~ 50%, since we include contributions to the 
amplitude up to next -to-leading order in (M^/A^) in our calculation. The factor of 2 appears 
since the width is proportional to the square of the amplitude. It is reassuring that both methods 
lead to essentially similar numbers for the uncertainty. Our final result therefore reads 

T (L>* (2317)+ -► D+tt°) = (180 ± 40 ± 100) keV, (24) 

where the first error is from experimental inputs and the second reflects the theoretical uncertainty. 
Amongst the former uncertainties, the largest ones are from the uncertainties of the .D-meson 
masses and the ir°-r] mixing parameter A = 0.024 ± 0.002. Schematically, let us consider the case 
corresponding to h' 5 = 1. The central result for the width is 233 keV. When the central values of 
all the meson masses are taken, the width can change from 219 keV for A = 0.022 to 248 keV for 
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300 




Figure 2: T (D* (2317) + — > -D+7T ) as a function of h' 5 . The solid curve results from choosing the 
central value for A, while the grey band is found by varying A within its allowed bounds and taking 
into account theoretical uncertainty. The dashed curve refers to neglecting the e.m. term in the 
decay (let g = -2g 2 ). 



A = 0.026. When we take the central value A = 0.024 and the central values of the masses of all 
the mesons except D° and D + , the width can change from 233 keV for taking the central values 
of m.£)0 and rri£,+ to 249 keV for taking m^o = 1864.65 MeV and m£)+ = 1869.82 MeV. Among 
all the others, the uncertainties caused by go + 2g 2 , see Eq. (16), and the masses of kaons are the 
largest, and they amount to an uncertainty of 3 keV and 1 keV at most, respectively. All other 
terms give negligible contributions to the uncertainty. 

Within the molecular picture, different calculations gave the width of the radiative decay 
D* (2317) — > D*"f in the range from 1 — 6 keV [27, 13, 14]. Combining with the experimental 
result in Eq. (1), the lower limit for the width of the decay -D* (2317) —* D s 7T° is of the order of 
100 keV. Our result is compatible with this extracted lower limit. 

5 Summary and Outlook 

In this paper we investigate the isospin violating decay D* (2317) — ► D s tt up to NLO in the 
chiral expansion, assuming that the -D* (2317) is a hadronic molecule. We take into account 
electromagnetic contributions systematically for the first time. Up to order 0(p 2 ), we obtain both 
the strong and e.m. mass differences of the Z)-mesons. We confirm that the mass differences 
between charged and neutral mesons in the same isospin multiplets play a significant role in the 
decay width. The decay width of the D* (2317) — > D s ir calculated to next-to-leading order is 
found to be 180 ± 110 keV, where the uncertainties are added in quadrature. The uncertainty is 
dominated by the theoretical one. The resulting width is consistent with the present experimental 
constraint. 

Our results for the hadronic decay width of the -D* (2317) within uncertainty are consistent with 
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previous analyses of the D* (2317) within a molecular picture considering both the Tr°-r] mixing 
and the meson mass differences: Ref. [13] gives 79.3 ± 32.6 keV, while Ref. [14] gives 76 keV (140 
keV) as the result at leading (next -to-leading) order. Assigning the D* (2317) to be a cs meson, its 
hadronic width was estimated within quark models with typical values of the order of 10 keV [28], 
although larger values were reported — see collection of results in table II of Ref. [13], which is 
consistent with the analysis utilizing heavy quark effective field theory [6] . Within the tetraquark 
picture a similar width as in the cs picture was found [29] , which is about one order of magnitude 
smaller than our predictions in the molecular picture. We therefore conclude that the decay width 
of the -D* (2317) — ► D s ir° can be a good criterion for testing the nature of the Z)* (2317). A 
simultaneous study of radiative decays within the various scenarios is also necessary, as advocated 
for the light scalar mesons in Ref. [30]. To expose the nature of the -D* (2317), experimental efforts 
are highly appreciated to improve the quantitative understanding of the D* (2317) decays. 

In this paper, we only considered the pseudoscalar D-mesons and the Goldstone bosons. The 
effects of all the higher states are incorporated in the LECs. However, the mass difference between 
D* and D is only about 140 MeV, which is approximately equal to (ttia — mi\i)/2. As it is 
sometimes important to include the A (1232) in the chiral effective field theory for baryons (for a 
recent review, see Ref. [18]), it would be interesting to check what would happen if we include the 
vector charm mesons explicitly in the effective Lagrangian, as in Ref. [14]. This extension of the 
scheme will be investigated in the future. 
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A Suppression of the Hq terms 

Due to the commutator structure in Eq. (6), all the h,Q terms in amplitudes are proportional to 

{Pi ■ P2)(P3 ■ Pa) ~ {Pi ■ Pa)(P2 ■ P3) , (25) 

where p\ (p 2 ) and p 3 (p^) are the momenta of the heavy mesons (Goldstone bosons) in the initial 
and final state, respectively. Let v denote the velocity of a heavy meson, we separate the momenta 
of the heavy mesons into two parts as 

Pi = Miv + ki, 

p 3 = M 3 v + k 3 , (26) 

where Mi and M 3 are the masses of the heavy mesons, and k\ and k 3 are small residual momenta 
which are of order O (p) . Let p 3 = Pi + Ap with Ap = (M 3 — Mi ) v + k 3 — k\ , we have 

{Pi -P2){P3 - Pa) ~ {Pi -Pa)(P2 -P3) = Ap - [{pi ■p 2 )pa ~ {Pi -Pa)P2] ■ (27) 

Because \M 3 — M\\ ~ 100 MeV at most, Ap should be counted as 0{p). Thus the above equation 
should be counted as 0(p 3 ), and hence is suppressed by one more order. 



9 



References 

[1] E. S. Swanson, Phys. Rept. 429 (2006) 243 [arXiv:hep-ph/0601110]; E. Klempt and A. Za- 
itsev, Phys. Rept. 454 (2007) 1 [arXiv:0708.4016 [hep-ph]]; M. B. Voloshin, arXiv:071 1.4556 
[hep-ph], Prog. Nucl. Part. Phys. in print. 

[2] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 90 (2003) 242001 [arXiv:hep- 
ex/0304021]. 

[3] W. M. Yao et al. [Particle Data Group], J. Phys. G 33 (2006) 1. 

[4] T. Barnes, F. E. Close and H. J. Lipkin, Phys. Rev. D 68 (2003) 054006 [arXiv:hep- 
ph/0305025]. 

[5] E. van Beveren and G. Rupp, Phys. Rev. Lett. 91 (2003) 012003 [arXiv:hep-ph/0305035]. 

[6] T. Mehen and R. P. Springer, Phys. Rev. D 72 (2005) 034006 [arXiv:hep-ph/0503134]; 
T. Mehen and R. P. Springer, Phys. Rev. D 70 (2004) 074014 [arXiv:hep-ph/0407181]. 

[7] C. Hanhart, Yu. S. Kalashnikova, A. E. Kudryavtsev and A. V. Nefediev, Phys. Rev. D 75 
(2007) 074015 [arXiv:hep-ph/0701214]. 

[8] F. K. Guo, S. Krewald and U.-G. MeiBner, arXiv:0712.2953 [hep-ph], Phys. Lett. B (2008) 
in print. 

[9] V. Baru, J. Haidenbauer, C. Hanhart, Yu. Kalashnikova and A. E. Kudryavtsev, Phys. Lett. 
B 586 (2004) 53 [arXiv:hep-ph/0308129]. 

[10] E. E. Kolomeitsev and M. F. M. Lutz, Phys. Lett. B 582 (2004) 39 [arXiv:hep-ph/0307133]; 
D. Gamermann, E. Oset, D. Strottman and M. J. Vicente Vacas, Phys. Rev. D 76 (2007) 
074016 [arXiv:hep-ph/0612179]. 

[11] J. Hofmann and M. F. M. Lutz, Nucl. Phys. A 733 (2004) 142 [arXiv:hep-ph/0308263]; 

[12] F. K. Guo, P. N. Shen, H. C. Chiang and R. G. Ping, Phys. Lett. B 641 (2006) 278 [arXiv:hep- 
ph/0603072]; F. K. Guo, P. N. Shen and H. C. Chiang, Phys. Lett. B 647 (2007) 133 
[arXiv:hep-ph/0610008] . 

[13] A. Faessler, T. Gutsche, V. E. Lyubovitskij and Y. L. Ma, Phys. Rev. D 76 (2007) 014005 
[arXiv:0705.0254 [hep-ph]]. 

[14] M. F. M. Lutz and M. Soyeur, arXiv:0710.1545 [hep-ph]. 

[15] S. Weinberg, Trans. New York Acad. Sci. 38 (1977) 185. 

[16] U.-G. Meifiner and S. Steininger, Phys. Lett. B 419 (1998) 403 [arXiv:hep-ph/9709453]; 
G. Miiller and U.-G. Meifiner, Nucl. Phys. B 556 (1999) 265 [arXiv:hep-ph/9903375]. 

[17] V. Bernard, N. Kaiser and U.-G. Meifiner, Int. J. Mod. Phys. E 4 (1995) 193 [arXiv:hep- 
ph/9501384]. 

[18] V. Bernard, Prog. Part. Nucl. Phys. 60 (2008) 82 [arXiv:0706.0312 [hep-ph]]. 

[19] G. Burdman and J. F. Donoghue, Phys. Lett. B 280 (1992) 287; M. B. Wise, Phys. Rev. D 
45 (1992) 2188; T. M. Yan, H. Y. Cheng, C. Y. Cheung, G. L. Lin, Y. C. Lin and H. L. Yu, 
Phys. Rev. D 46 (1992) 1148 [Erratum-ibid. D 55 (1997) 5851]. 



10 



[20] H. Y. Cheng, C. Y. Cheung, G. L. Lin, Y. C. Lin, T. M. Yan and H. L. Yu, Phys. Rev. D 
49 (1994) 5857 [Erratum-ibid. D 55 (1997) 5851] [arXiv:hep-ph/9312304]. 

[21] A. V. Manohar, arXiv:hep-ph/9802419. 

[22] A. Kucukarslan and U.-G. Meifiner, Mod. Phys. Lett. A 21 (2006) 1423 [arXiv:hep- 
ph/0603061]. 

[23] N. Fettes and U.-G. Meifiner, Phys. Rev. C 63 (2001) 045201 [arXiv:hep-ph/0008181]. 
[24] J. Gasser and H. Leutwyler, Phys. Rept. 87 (1982) 77. 

[25] J. A. Oiler and U.-G. Meifiner, Phys. Lett. B 500 (2001) 263 [arXiv:hep-ph/0011146]. 

[26] J. A. Oiler and E. Oset, Phys. Rev. D 60 (1999) 074023 [arXiv:hep-ph/9809337]. 

[27] D. Gamermann, L. R. Dai and E. Oset, Phys. Rev. C 76 (2007) 055205 [arXiv:0709.2339 
[hep-ph]]. 

[28] S. Godfrey, Phys. Lett. B 568 (2003) 254 [arXiv:hep-ph/0305122]. 
[29] M. Nielsen, Phys. Lett. B 634 (2006) 35 [arXiv:hep-ph/05 10277]. 

[30] Yu. Kalashnikova, A. E. Kudryavtsev, A. V. Nefediev, J. Haidenbauer and C. Hanhart, Phys. 
Rev. C 73 (2006) 045203 [arXiv:nucl-th/0512028]. 



11 



